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Next-to-leading order QCD corrections to ttZ production at the LHC 
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We present a calculation of the full next-to-leading order QCD corrections to the scattering 
process pp —> ttZ. This channel will be used to measure the ttZ electroweak couplings at the Large 
Hadron Collider. These couplings cannot be directly measured in current experiments. To obtain 
these results, we utilize a novel, completely numerical approach to performing higher order QCD 
calculations. We find that for reasonable values of the renormalization and factorization scales 
the QCD corrections increase the leading order result by 35%, while reducing the theoretical error 
arising from uncalculated higher order corrections to approximately ±10%. 



Probing the properties of the top quark is an urgent 
priority of the particle physics experimental program. 
The large mass of the top quark, more than an order of 
magnitude larger than any other fermion, suggests that it 
is intimately connected to the mechanism of electroweak 
symmetry breaking and to the generation of hierarchies 
in the flavor sector of the Standard Model (SM). Al- 
though it was discovered more than a decade ago at the 
Fermilab Tevatron 1] , relatively little is known about the 
top quark. The production process pp — > tt at the Teva- 
tron proceeds almost entirely through giuon exchange, 
and is sensitive only to the top quark mass and SU {i)c 
representation. The top quark decays primarily through 
t — » Wh, leaving only the helicity of the final state W to 
be measured. 

Details regarding the electroweak properties of the top 
quark are unknown. This information requires measure- 
ment of the ^7, ttZ , and top quark Yukawa couplings. 
These couplings probe many new physics effects; they are 
sensitive to mixing with additional Z' gauge bosons and 
new heavy fermions, and can provide access to extended 
scalar sectors. Although inaccessible at the Tevatron, 
these interactions can be studied through the radiative 
processes pp — > ti^^ Z, H at the Large Hadron Collider 
(LHC). Measuring deviations from SM predictions can 
provide further information regarding new states discov- 
ered at the LHC, or can provide evidence of new physics 
too heavy to produce directly. 

The cross sections for the scattering channels pp 
ttj, Z, H scale as (Tlo ^ o?s leading order in the QCD 
perturbative expansion, indicating that there is signif- 
icant theoretical uncertainty in predicting their rates. 
It is important to quantify the effect of this error on 
the measurement of the top quark couplings and deter- 
mine whether a next-to-leading order QCD calculation 
is needed to improve the theoretical prediction. A pro- 
cedure for measuring the tfy, Z couplings was discussed 
in 0. This study showed that the uncertainties in the 
signal cross sections hinder the extraction of anomalous 
top quark couplings. The theoretical error coming from 



higher order QCD corrections is the limiting factor in 
measuring ttZ couplings at the LHC. 

In this Letter we present a full calculation of the next- 
to-leading order (NLO) QCD corrections to the pp ttZ 
scattering channel to facilitate accurate measurements of 
top quark properties at the LHC. Computations of NLO 
QCD corrections to processes with five or more exter- 
nal particles are notoriously difficult, and require over- 
coming many technical challenges. Significant commu- 
nity effort has been invested in devising efficient calcula- 
tional algorithms for one- loop calculations ^ , leading to 
a host of new results of phenomenological relevance for 
the LHC 13 . 

Recently, a completely numerical and highly flexible 
approach to NLO and NNLO calculations was suggested 
[all) 0, Hi • It applies sector decomposition @ and con- 
tour deformation [10| to Feynman parametric loop in- 
tegral representations. It avoids many pitfalls present in 
traditional approaches to calculations of QCD corrections 
to multi-particle production processes. The calculation 
of the QCD corrections to pp ttZ described here is the 
first time these techniques have been applied to a com- 
plete scattering channel with multiple mass scales. This 
mode contains the full spectrum of difficulties present 
in the most complex 2 — > 3 processes. Several five-point 
topologies with two, three and four internal massive prop- 
agators appear. In this Letter we outline the technical 
details of our calculation and present phenomenological 
results relevant for LHC analyses. Preliminary results for 
the gg initiated partonic channel were presented in [6j. 

We consider the process p{Pi) + p{P2) ttZ. The 
factorization theorems for hard scattering processes in 
QCD allow us to write 

1 

dcr==^y" dxidx2flixi)ff{x2)daij^ttz{xiX2S), (1) 

f 1 2) 

where f-j (2^1,2) are the parton densities that give the 
probability to find parton in the proton 1(2) with 
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momentum = 2^1(2) ^'1(2)- The center-of-mass energy 
squared of the proton-proton coUision is introduced in 

Eq. d), 5 = 2Pi - Pa. 

At leading order in the as expansion, both gg and qq 
initiated partonic processes contribute. The computa- 
tion of the leading order cross section is straightforward. 
We use QGRAF [11] to generate the relevant Feynman 
diagrams and then MAPLE and FORM [l^ to manip- 
ulate this output. Throughout the paper we set nit = 
170.9 GeV, mz = 91.19 GeV, and mw = 80.45 GeV. 
For the coupling of the Z boson to quarks, we employ 

Zqq: iJ ^f^^/ {gl + gll,) , (2) 
y V 2 cos^ Ow 

rpq rpq 

where gl ^ - Qq sir? Ow, 5a = --|-, sin^ Ow = I - 

m^/m?^ — 0.2215 is the sine squared of the electroweak 
mixing angle, Tg is the weak isospin of the quark q, Qq is 
the electric charge of the quark q in units of the proton 
charge and Gi? is the Fermi constant. Numerical results 
for the leading order cross section are presented below. 

At next-to-leading order, the qg channel additionally 
contributes. Several distinct contributions to the NLO 
cross section can be identified: 

• the one-loop virtual corrections to the leading order 
processes gg, qq ttZ; 

• the real emission corrections gg — > ttZg, qq — > ttZg 
and qg — *■ ttZq\ 

• the renormalization of the leading order cross sec- 
tion. 

Both the virtual and real emission corrections are sep- 
arately divergent. They must be combined to remove 
divergences arising from soft gluon emission. Ultraviolet 
divergences are absorbed into the definitions of the cou- 
pling constant, the top quark mass, and the top quark 
wavefunction. Singularities associated with initial state 
coUinear gluon emission are absorbed into the definition 
of the parton distribution functions. We renormalize 
the top quark mass and wavefunction on shell. We em- 
ploy the MS scheme for the parton distribution functions, 
and for the QGD coupling constant as for the five light 
quarks. The ultraviolet divergence coming from the top 
quark is removed through a zero-momentum subtraction. 
These renormalizations are performed at certain momen- 
tum scales. If computed to all orders in perturbation the- 
ory the cross section would be independent of this scale 
choice; the residual dependence of the NLO result on the 
scale provides an estimate of the uncalculated higher or- 
der corrections. 

We utilize dimensional regularization to control the 
singularities of the virtual and radiative corrections at 
intermediate stages. The radiative corrections are orga- 
nized using the two-cutoff slicing method [Tsj . To com- 
pute the one-loop virtual corrections we employ the strat- 
egy we introduced in 0, @1 ■ For each one- loop diagram 



interfered with the full Born amplitude, we analytically 
perform the integration over the loop momentum and 
arrive at a Feynman parametric representation. While 
the ultraviolet divergences factorize after the integration 
over the loop momentum and therefore can be simply 
extracted, the infrared and coUinear singularities appear 
on the boundaries of Feynman-parameter space. We 
must therefore extract the infrared and coUinear singular- 
ities before integrating numerically. To do so we sector- 
decompose the Feynman-parametric integrals. However, 
even after these divergences are extracted, it is not possi- 
ble to perform the Feynman-parametric integrations nu- 
merically because there are singularities inside the in- 
tegration region coming from internal loop thresholds. 
These singularities are avoided by deforming the inte- 
gration contour into the complex plane. Once infrared 
and coUinear singularities are extracted and the integra- 
tion contour is deformed, we obtain representations of the 
Feynman parametric integrals suitable for numerical in- 
tegration. No reduction of tensor integrals is performed. 
This allows us to bypass many of the issues associated 
when this algebraic reduction is attempted. 

We present below phenomenological results of our 
NLO QGD calculation of ttZ . We employ 

the Martin-Roberts-Stirling-Thorne parton distribution 
functions [Ol at the appropriate order in the perturba- 
tive expansion. All numerical results are obtained us- 
ing the adaptive Monte Garlo integration algorithm VE- 
GAS [Hi as implemented in the GUBA library [li|. We 
have checked that our results satisfy several consistency 
checks. The leading order result obtained with our code 
matches that obtained using the program MadEvent [13] . 
Using an eikonal approximation, the divergent parts of 
the one-loop virtual corrections arising from soft gluon 
exchange can be calculated in a simple analytic form. 
We have checked that our full results agree with this 
form. All divergences cancel once we assemble the sep- 
arate components of the computation discussed above. 
The virtual corrections are independent of the size of 
our contour deformation of the Feynman parametric in- 
tegrals. Finally, we have implemented all parts of our 
calculation in several independent codes which agree for 
all observables studied. 

We present in Fig. ^ the inclusive cross section at 
both leading order and next-to-leading order in the per- 
turbative expansion. We have equated the renormaliza- 
tion and factorization scales to a common value /i — 
M-R — /^F, and have varied them from /io/8 to 2/io with 
/iQ — 2'mt + mz. The dependence of the pp ttZ rate on 
this unphysical scale parameter is significantly lessened 
when the NLO corrections are included. Choosing a scale 
too different from the typical momenta and energies that 
give the dominant contributions to the process leads to 
large logarithms that spoil the convergence of the pertur- 
bative expansion. Therefore, to estimate a central value 
and theoretical error for the cross section, we should a 
pick a scale roughly equal to the typical transverse mo- 
menta and masses in the final state and vary around 
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this value. As these momenta and masses are approxi- 
mately 100-200 GeV, we consider in [/xo/4,/io] a rea- 
sonable range of scale variation with /i = /io/2 a good 
central value. This yields a cross section of 1.09 pb with 
a theoretical error of ±11% at NLO. The result at LO is 
0.808 pb with an uncertainty of ±25 — 35%. The inclu- 
sive Kinc-fstCtoT for this process, defined as the ratio of 
the cross section at NLO to that at LO, is Kinc = 1.35 for 
fi ~ Mo/2. The variation of Kmc with scale is also shown 
in Fig. Il]); it changes from 1.1 to 1.6 as /i varies from 
/j,o/4 to /iQ. Also included in Fig. ([1]) are the separate 
contributions of the gg, qg, and qq partonic processes at 
NLO. The significant scale dependence of the qg compo- 
nent, which first appears at this order in the perturbative 
expansion, is noteworthy. 



the tt pair. This intuitive argument leads us to expect 
that the shape of many other kinematic distributions will 
also be approximately unchanged by NLO corrections. 
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FIG. 1: Inclusive cross section for pp ttZ as a function of 
the scale choice fi. Included are the LO and NLO results, as 
well as the contributions of the gg, qg, and qg partonic chan- 
nels at NLO. The dotted line shows the inclusive A"- factor; 
the value for this should be read from the axis on the right of 
the plot. 

In addition to the inclusive cross section, the impact 
of higher order corrections on differential distributions 
must be studied. An interesting question to consider is 
whether their effect is completely described by the in- 
clusive Xinc-factor. If so, the NLO corrections can be 
accurately and simply included in leading order simula- 
tion codes by an overall reweighting of event rates. To 
investigate this question we present in Fig. ^ the bin- 
integrated transverse momentum spectrum of the Z bo- 
son at both LO and NLO for the scale choice /i = fJ,o/2. 
Most of the cross section comes from events with less 
than 200 GeV. Included in this plot is the ratio of the 
NLO distribution over the LO spectrum, Kp^. It is 
fiat to within a few percent over the entire range, and 
is equal to the inclusive value Kinc — 1-35. The small 
impact of higher order corrections on the distribution 
can be roughly understood by noting that at tree level, 
pp — > ttZ is already a three-body process. Including ad- 
ditional partonic radiation does not open up new regions 
of phase space as the Z boson can already recoil against 
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FIG. 2: Transverse momentum spectrum for pp tiZ for the 
scale choice /i = /io/2 = mt + mz/'2. Included are the LO and 
NLO results. The pt dependent if -factor for each bin, Kp^, 
is also shown; the value for this should be read from the axis 
on the right of the plot. 

We can use these results to estimate the improvement 
in the measurement of tiZ couplings at the LHC af- 
ter NLO corrections are included. Assuming approxi- 
mate CP conservation, four relevant tiZ couplings ex- 
ist: dimension-four vector and axial couplings, and two 
dimension-five dipole couplings. Although the measure- 
ments of these parameters are correlated, the analysis of 
Ref. f3| indicates that the dipole couplings are expected 
to be measured with a precision of ±50% at the LHC 
and ±25% with the super-LHC luminosity upgrade, while 
the axial coupling should be measured with ±15% preci- 
sion at the LHC and with ±5% at the super-LHC. This 
study also found that the uncertainty arises almost en- 
tirely from the signal normalization and statistics; the 
backgrounds are negligible. This analysis utilized a scale 
choice /i = mt, for which we find Kinc ~ 1.3. This yields 
a small improvement of the relative statistical error. The 
theoretical uncertainty assumed in this analysis of Q was 
±30%. The authors also studied the expected improve- 
ment possible if higher order corrections reduced this er- 
ror to ±10%, and concluded that improvements in the 
precisions quoted above could reach a factor of two at the 
LHC. A conservative estimate of the remaining theory 
uncertainty from our prediction for pp — > ttZ is ±15%. 
This accounts for imprecise knowledge of parton distribu- 
tion functions. While the full factor-of-two improvement 
in the precision of the ttZ couplings seems slightly out 
of reach with our current knowledge of the higher order 
corrections, we still expect a significant improvment once 
this reduced theoretical error is propagated through the 
full analysis. Presumably the improvement is more sig- 
nificant with the super-LHC luminosity upgrade, since 
the relative importance of the statistical errors should be 
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decreased. 

In summary, we have presented the complete calcula- 
tion of the next-to-leading order QCD corrections to the 
pp ttZ scattering process at the LHC. Study of this 
channel allows measurements of the ttZ electroweak cou- 
plings, which cannot be directly accessed with current ex- 
periments. These couplings probe many forms of physics 
beyond the Standard Model, such as small mixings of ex- 
tra heavy gauge bosons or vector-like fermions with the 
Standard Model top quark and Z boson. To compute 
these corrections we present a novel approach to pertur- 
bative calculations in QCD powerful enough to handle 
the significant complexities that occur when 2 — > 3 and 
more complicated scattering processes are studied. This 
method is highly automated and flexible, and is based 
on a completely numerical algorithm for computing loop 
integrals. This is the first time it has been applied to a 
scattering process that contains the full set of complexi- 
ties that occur in multi-leg one-loop calculations. We find 
that the NLO QCD corrections increase the leading order 
pp ttZ cross section by a factor of 1.35 for canonical 
choices of the factorization and renormalization scales. 
We estimate the remaining theoretical uncertainty from 



uncalculated higher order corrections to be ±11%. Pre- 
vious studies based on leading order cross sections have 
found that the normalization uncertainty of the pp — > ttZ 
rate is the dominant error in extracting ttZ couplings at 
the LHC. We estimate that the relative precision of the 
ttZ couplings will be improved by a factor of 1.5-2 once 
our results are incorporated into these analyses. We find 
that the NLO corrections do not significantly change the 
shape of the kinematic distributions we have studied, in- 
dicating that they can be accounted for by an overall 
scaling of the LO result. 

We are excited by the potential of our approach to 
NLO QCD calculations presented here. We anticipate 
and look forward to using it to understand other pro- 
cesses of interest at future collider experiments. 
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